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UNCERTAINTY IN DETERMINING THERMAL PERFORMANCE OF LIQUID-HEATING FLAT-PLATE SOLAR COLLECTORS 
By 


E.R. Streed and D. Waksman 


ABSTRACT 


Thermal performance measurements of eight liquid-heating flat-plate solar collectors were 
conducted with two to four collectors of each type at four outdoor test sites. Tests were 
performed in accordance with the procedure prescribed by ASHRAE Standard 93-77. Statistical 
analysis of data sets for each collector type within test sites and between test sites was 
done using ASTM recommended methods to evaluate test method measurement uncertainty. 
Illustrations of the influence of thermal performance data uncertainty are presented for 
collector material degradation, collector rating and calculated system performance. 


KEY WORDS: Collector rating; measurement; solar collector; standards; thermal 
performance; uncertainty. 
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PREFACE 


The Department of Energy (DoE) is conducting an Interim Solar Collector Testing Program to 
stimulate the establishment of a permanent collector testing/certification program. 

Elements of the program are related to\the verification of test procedures, evaluation and 
improvement of test laboratory capabilities, and the measurement and distribution of . 
reliable performance data for solar collectors. This study of collector thermal performance 
data uncertainty was conducted by the National Bureau of Standards (NBS) in support of this 
DoE test program. 


This report is intended to be informative and instructive; it is not meant to be an 
evaluation of products or facilities. In no case does such identification imply 
recommendation or endorsement by the National Bureau of Standards, nor does it imply 
that the products are necessarily the best available for the purpose. 


The cooperation of the following people is greatly appreciated: Richard D. Whitaker and 
William Dokos, DSET Laboratories, Inc.; Ross McCluny, James Huggins and James Rowland, 
Florida Solar Energy Center; Roger Wedell and Ronald Dammann, Lockheed Palo Alto Research 
Laboratory; Ken Shih and Mark Rallas, Wyle Laboratories; and Aaron Dawson, Virginia 
Polytechnic Institute and State University. The devoted effort of Dr. John Mandel in per- 
forming the statistical analysis, instrument calibration by John Jenkins, data collection 
by Cathy Scarbrough and Patricia Christopher, and computerized calculations by William May, 
all at NBS, contributed significantly to the evaluation and presentation of the data. 
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SI CONVERSION UNITS 


In view of the present accepted practice in this country for building technology, common 
U.S. units of measurement have been used in the data appendix. In recognition of the 
position of the United States as a signatory to the General Conference of Weights and 
Measures, which gave official status to the metric SI system of units in 1960, assistance 

is given to the reader interested in making use of the coherent system of SI units by giving 
conversion factors applicable to U.S. units used in this document. 


Length 
1 in = 0.0254 meter (exactly) 


1 ft = 0.3048 meter (exactly) 
Area 

1 in* = 6.45 x 107" meter? 

1 £t? = 0.09290 meter? 
Volume 


1 in = 1.639 x 107° meter? 
1 gal (U.S. liquid) = 3.785 x 10-2 meter? 


Mass 
1 ounce-mass (avoirdupois) 2.834 x 10-* kilogram 
1 pound-mass (avoirdupois) = 0.4536 kilogram 


Pressure or Stress (Force/Area) 
1 inch of mercury (60°F) = 3.377 x 102 pascal 
1 pound-force/inch* (psi) = 6.895 x 10° pascal 


Energy 
1 foot-pound-force (ft-lbf) = 1.356 joule 
1 Btu (International Table) = 1.055 x 102 joule 
1 kilowatt-hour = 3.600 x 10© joule = 3.412 x 10° Btu 


Power 
1 watt = 1 x 10’ erg/second 
1 Btu/hr = 0.2929 watt 


Temperature 
foo = 5/9 (top ae 32) 


Heat 
1 Btu-in/h:ft2-°F = 1.442 x 10-! W/m-K (thermal conductivity) 
1 Btu/lbm-°F = 4.184 x 10° J/kg-K (specific heat) 
1 langley = 4.184 x 10% J/m2 = 1 cal/cm* = 3.69 Btu/ft2 
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INTRODUCTION 
1.1 Background 


The Department of Energy (DoE) and the Department of Housing and Urban Development 
(HUD) are committed to the demonstration of building solar heating and cooling 
applications in the United States. One of the objectives of the demonstration 
program involves the development of standardized procedures and methods for testing 
and evaluating solar collectors. As stated in the "National Program for Solar 
Heating and Cooling of Buildings" [1], the National Bureau of Standards (NBS) is 
working with industry and professional organizations to develop nationally recog- 
nized consensus test procedures. 


The ASHRAE Standard 93-77 [2] has been adopted in the United States (U.S.) as a 
consensus test method for determining the thermal performance of solar collectors 
and is generally recognized as an adequate procedure. Initial development of the 
test procedure [3] and experimental evaluation for selected flat-plate, evacuated 
tubular and low concentrating collectors [4] have been reported by NBS. A round 
robin type study of the basic test procedure was conducted utilizing 21 test 
facilities with two types of flat-plate collectors during 1976 and 1977. Analysis 
of these results [5] showed that correction for the influence of the environment 
and test conditions to a set of reference conditions for incidence angle, diffuse 
fraction, wind, irradiance and ambient temperature resulted in only a minor reduc- 
tion in data spread. It was concluded that systematic type measurement problems 
resulting from instrument calibration, sensor installation and procedural varia- 
tions are probably more significant contributors to the data uncertainty. 


1.2 Current Studies 


In order to stimulate the establishment of an industry operated collector testing and 


certification program, DoE organized an Interim Collector Testing Program in 1978. 
The program had three major objectives as follows: 


1. Provide state-of-the-art thermal performance data for designer, 
manufacturer and consumer use on a large number of collectors 

2. Provide a technical base for the establishment of a collector 
certification and rating program 

3. Evaluate and modify (if necessary) collector performance test 
methods 


The thermal performance measurements portion of the program has been completed and 
the data are being compiled for publication. Only a limited number of tests were 
conducted for comparability purposes. 


The primary factors that need to be determined are a measure of the repeatability 
and reproducibility of ASHRAE Standard 93-77. Repeatability is related to the 
variability in replicate results obtained on the same product within a single test 
facility. Reproducibility is related to the variability in results obtained on 
the same product between test sites. 


As part of the overall DoE Solar Energy Program, a project was initiated in 1978 by 
NBS [6] to investigate accelerated and "real time'’ test methods for predicting 
collector durability and reliability. Representative commercially available flat- 
plate solar collectors and their materials are being exposed outdoors in four 
different U.S. climatic regions and with two types of solar simulators. As part 

of the characterization of the collectors, initial thermal performance measurements 
were made on multiple samples of each collector at each test site prior to 
commencing the exposure tests. These baseline data consisting of measurements 

made on eight collector types at each of four sites, with two to four collector 
replicates at each site, provide an appropriate data set which is used in this 
study to evaluate collector thermal performance measurement capabilities. 


The NBS test programwas primarily designed for durability/reliability considerations, 
and does not include some features desirable for identifying individual factors 
contributing to the measurement uncertainty, such as product variability, specific 
environmental conditions, and test apparatus. Additional investigations to evaluate 
and improve the ASHRAE test method are currently being conducted under the direction 
of an ASHRAE 93R (Revision) Committee. 


1.3 Report Purpose 


It is the intent of this study to show the following: (1) overall uncertainty 

in collector thermal performance data as reported by commercial test facilities 

on commercially available solar collectors (individual error sources will not be 
identified except for some cursory observations of possible causes), and 

(2) illustrations of how’the uncertainty in individual collector test results affects 
the rating and selection of collectors for system design and for determining degrada- 
tion with exposure or operating time. 


THERMAL PERFORMANCE TEST METHOD 
2.1 General 


The thermal performance tests used to obtain the instantaneous efficiency were 
specified to be performed in general accordance with ASHRAE Standard 93-77. This 
consensus test method for both air and liquid heat transfer fluids includes guide- 
lines for the test apparatus configuration, instrumentation accuracy, and precision 
tolerances. The measurements needed to determine the instantaneous efficiency are 
temperature, fluid mass flow, solar irradiance, and collector area,as expressed by 
the following equation; 


n= G. (core - eee (1) 

where fn - = collector efficiency 

m = fluid mass flow rate (raison) 

c = fluid specific heat (J/kg: °C) 

cet = transfer fluid inlet temperature (°C) 

Ceo = transfer fluid exit Sa meani GG) 

A = collector gross area (m ) 

i. = solar irradiance Gime. 


These measurements must be performed precisely with calibrated instruments capable of 
achieving the prescribed accuracy and also must be performed within a limited range 
of weather ana test conditions. During testing, the solar irradiance, inlet 
temperature, and fluid flow rate must be maintained within close tolerances for time 
periods in excess of the collector time constant to insure that a quasi-steady state 
thermodynamic condition prevails in the collector with respect to the gain and loss 
of energy. Climatic factors such as ambient temperature and wind have relatively 
wide allowable ranges and they must be measured and reported as well as the solar 
diffuse fraction and relative humidity. Other test parameters which may influence 
the thermal performance but are not specifically restricted are tilt angle, sky 
temperature and solar spectral distribution. 


The test procedure specifies that at least four inlet temperatures shall be used to 
obtain near normal incidence efficiency values over a range of operating conditions. 
At least four data points shall be taken at each inlet temperature. When measurements 
are made outdoors on a non-tracking mount, two data points should be taken before 
solar noon and two after solar noon to cancel out heat capacity effects. Measurements 
made using a solar simulator or an altazimuth mount are conducted with the sun 

normal to the collector during the test period. Each data point is obtained by 
integrating the measured parameters over a five-minute period (longer for collectors 
with time constants in excess of five minutes) and calculating the efficiency as 
indicated in equation 1. 


2.2 Efficiency 
ASHRAE Standard 93-77 is based upon the Hottel-Whillier-Bliss analytical model which 


assumes that the efficiency of a flat-plate solar collector operating under steady- 
state conditions and normal incidence can be described by: 


n= F(t) - FRU (te 5 = t/t] (2) 


where FR = collector heat removal factor 
(ta) , = effective transmittance-absorptance product 
UL = heat transfer loss coefficient [W/(m2°C) ] ; 


The concept of the collector heat removal factor, Fp, was introduced in the testing 
of thermal collectors to preclude the necessity for determining the mean temperature 
of the absorber plate. The parameter Fp is defined as: 


rae actual useful energy collected by a flat-plate collector 
R useful energy collected if the entire flat-plate collector 


surface were at the inlet fluid temperature 
The quantity Fp may be expressed in exponential form as follows: 
th c 
= F! 
% CU, A,/m ey 


= = (3) 
ROA 


F 


2 
where A, = collector aperture area (m). 


The collector efficiency factor F"' is made up of a series of heat path resistances 
dependent upon the absorber physical construction (materials and configuration) and 
the heat transfer coefficient from absorber surface-to-heat transfer fluid. Thus, 
the effects of collector physical design, mass flow rate, and fluid thermophysical 
properties are manifested in the collector's measured efficiency, assuming that the 
external boundary conditions of total incident radiation flux Gas ambient 
temperature (t,)> sky temperature, wind velocity, etc., are relatively constant. 


A plot of the collector efficiency versus the AT/I; parameter for the range of 
inlet temperatures can result in a linear relationship. The intercept on the "y" 
axis is equal to Aa/Ag Fr(ta)~* and the slope js equal to Aa/Ag FRU;*. Because 

of non-linearity in the loss coefficient with temperature, U;, is not always con- 
stant, and the test method provides for presenting the data by either a first-order 
or a second-order polynomial fit to the data points. 


A typical plot of the efficiency curves for a collector with a double cover and a 
flat black absorber and for one with a single cover and a selective absorber are 
shown in Figure 1. Examples of variation in wind and diffuse fraction are shown 
for individual data points obtained on different measurement days. The influence 
of environmental conditions on individual data points is smoothed out by fitting a 
curve to all 16 of the measured data points. 


2.3 Incident Angle Modifier 


ASHRAE Standard 93-77 also provides for measurements of the influence of solar 
incidence angle on the collector efficiency. The incident angle modifier, Kg-,is 
used in this method as a coefficient of the transmittance-absorptance product, 

(ca) 3 = Kr (ta), mp =o describe non-tracking collector performance when irradiated 
off-normal. The coefficient is used in calculating collector all-day performance. 
The value of Ky; for many flat-plate collectors can be plotted as a function of 


*The reported values are based upon calculations using collector gross area. Each term 
should be multiplied by the ratio of the gross area-to-aperture area (Ag/Aa) to conform 


to the definition of the heat removal factor expressed in equation 3, which is based upon | 


aperture area. 
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(1/cos 9)-1 as shown in Figure 2 and typically yields a straight line. The values 


of Kyr are determined by performing separate efficiency tests with the inlet 
temperature controlled to within +1°C (+2 °F) of the ambient temperature and at 
incident angles of 0, 30, 45 and 60°. This linear relationship shown in Figure 2 
may not hold for collectors with special (non-flat) geometry, reflectors, etc. 
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THERMAL PERFORMANCE DATA UNCERTAINTY 
3.1 General 


The evaluation of a test method involves an analysis of the total measurement 
process including: 


measurement personnel expertise 
test facility quality 

local environmental factors 
test specimen characteristics. 


0000 


The identification of individual sources of uncertainties would require the use of 
well characterized reference samples for measurement under identical environmental 
conditions or the exchange of the same specimen for measurements by each test 
facility. The conduct of this durability program did not permit sample control to 
this extent. The evaluation was conducted in a manner similar to the procedures 
employed by industry in a certification program. Collectors are typically randomly 
selected from production lots, sent to a qualified laboratory and tested under the 
prevailing environmental conditions within the prescribed test limitations. 
Therefore, typical variations in collector production, pre-test collector handling, 
test apparatus calibration, test procedure,and data reduction are aggregated in the 
data. 


Some factors which contribute to the uncertainty can be determined from analysis 

of the information obtained as part of the test procedure. For example, the 
collector gross area must be measured for use in calculating the efficiency as 
expressed by equation 1. Analysis of the reported area data for Test Site No. 1 
indicates a coefficient of variation cv; of 0.33 percent for "within" a test site. 
Similar analysis for the four test sites indicates a coefficient of variation cvR 
of 0.94 percent for "between" test sites. Assuming that the "within" value reflects 
the average variation in manufacturing tolerances for these collectors and that the 
larger "between" value is more indicative of measurement uncertainties, then the 
area determination becomes a significant part of the coefficient of variation for 
intercepts calculated in Section 5.2. 


The area, standard deviation, and coefficient of variation data are tabulated in 
Tables B-2 and B-3 of Appendix B. 


3.2 Measurement Uncertainties 


The influence of individual and cumulative instrumentation calibration and opera- 
tional errors has been calculated [7] using propagated measurement uncertainty 
functions. Uncertainty functions were determined for the allowable accuracies 
specified in the test procedure (assuming 2—o instrumentation errors) for 
individual measurements of solar irradiance, fluid flow rate, temperature and working 
fluid thermophysical properties plus a 1 percent uncertainty associated with the 
data acquisition. The magnitude of the caiculated measurement uncertainty asso- 
ciated with parameters comprising the x and y axes is illustrated in Figure 3 by 
the rectangles superimposed on the solid lines. Uncertainties of about +4 
percentage points are indicated for typical flat-plate collectors. The largest 
assumed individual contributing parameter is the uncertainty of +3 percent related 
to the measurement of solar irradiance. 
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3.3 Meteorological Uncertainty 


Additional variability in efficiency will result from the meteorological variations 
allowed in the test method. The range of allowable test conditions for solar 
irradiance, incidence angle, ambient temperature and wind speed specified in the test 
method are shown in Table 1. 


Table 1. Environmental Test Conditions Allowed by ASHRAE Standard 93-77 


Environmental Parameter 


Ambient air temperature range < 30°C 

Wind speed across collector should be < 4.5 m/s 
Total solar irradiance within collector plane > 630 W/m2 

Beam solar irradiance incident angle <0 

Reflected solar irradiation < 0.20 


The difference in calculated thermal performance for a single glass cover, 

selective absorber type collector exposed to these extreme conditions is represented 
by the solid curves shown in Figure 3. Calculations were made using an analytical 
model [5] based upon traditional heat transfer relationships for flat-plate 
collectors. It should be noted that the simultaneous occurrence of these extreme 
"high" or "low" meteorological conditions is extremely unlikely. 


Experimental verification of the calculated meteorological extremes was performed 

in a solar simulator to the extent that the conditions could be achieved. The solid 
black squares represent data points taken using the extreme exposure conditions 
listed under "simulator tests'' in Figure 3. The experimental data points are 
bracketed around the solid curve labeled "average outdoor}' which represents the mean 
value reported for round-robin type tests [5] on a similar collector. The experi- 
mental data were obtained with a xenon type solar simulator [8]. 
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BASELINE DATA 


4.1 Collector Description 


Liquid-heating flat-plate type collectors were selected for use in the test program. 
The designs chosen were representative of commonly used materials and types of 
construction. All collectors from each manufacturer were from the same production 
lot. The cover and absorber materials, their pertinent optical properties and 
average areas (gross and aperture) are shown for the respective collectors in 

Table 2. The absorber material optical property data are based upon measurements 

of at least ten samples taken from the actual absorber of each collector type. The 
solar transmittance of the glass materials was obtained using the full cover and a 
pyranometer (ASTM Standard E424 Method B [9]). The solar transmittance of non-glass 
cover materials and the solar absorptance values were obtained from spectral measure- 
ments with an integrating sphere (ASTM Standard E424 Method A [9]). Emittance was 
measured using a portable-type instrument employing a thermopile and infrared 
reflectance technique,in accordance with ASTM Standard E408, Method A [10]. 


Detailed descriptions of each collector's construction dimensions and pertinent 
material properties useful for thermal analytical modeling are listed in Table B-l, 
Appendix B. 


4.2 Test Sites 


One of the major objectives of the durability program is to determine the influence 
of environmental exposure conditions on collector degradation. The same type of 
information is required for evaluation of thermal performance testing. 


Four exposure test sites were selected which represent both the median and extreme 
United States climatological conditions. These test sites can be briefly described 
as follows: 


Climatological Extremes 


Sate hot. diy; 
high solar radiation 
(high UV radiation would accompany these conditions) 


The hot, dry condition can be found in southwestern states 
(i.e. Arizona, Nevada and New Mexico). DSET Laboratories, Inc., 
located in Phoenix, Arizona was selected as this test site. 


Site 2. hot, humid 
high solar radiation 
(low to moderate UV radiation would accompany these conditions) 


The hot, humid condition can be found either in Florida or along the 
Gulf Coast in the states of Alabama, Mississippi, Louisiana or Texas. 
The Florida Solar Energy Center, located in Cape Canaveral, Florida was 
selected as this test site. 


Median Climatological Conditions 


Site 3. moderate temperature 
dry 
high solar radiation 
moderate UV radiation 


The moderate, dry condition can be found primarily in California. The 


Lockheed Research Laboratory, located in Palo Alto, California, was 
selected as this test site. 
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Site 4. moderate temperature 
humid 
moderate solar radiation 
moderate to low UV radiation 


The moderate, humid condition can be found in the Pacific Northwest, 
Mid-Atlantic and Mid-South regions of the United States. The NBS test 
facility, located in Gaithersburg, Maryland will serve as this test site. 


Each of the commercial sites was determined by an industry inspection team to be 
capable of conducting tests in accordance with ASHRAE Standard 93-77 [11]. In 
addition, on-site checks of the instrumentation used by these laboratories were 
made by the Solar Energy Research and Education Foundation (SEREF) with their 
Laboratory Intercorrelation Test Equipment (LITE) apparatus as part of the DoE 
Collector Testing Program. In general, they were found to meet the required 
ASHRAE measurement specifications. Since the main concern of this study was in 
establishing the precision of the test method, detailed study of the type and 
accuracy of the individual sensors and apparatus used at each test facility was 
not considered necessary. 


4.3 Measured Data 


Four test series, each having a sample of the eight collector types, were measured 

at Sites 1 and 2,and two test series at Sites 3 and 4. The thermal performance of 
all collectors in Series 1, 3 and 4 was measured prior to outdoor stagnation 

testing. Collectors in Series 2 were exposed for three days to a minimum daily 

solar radiation of 17,100 KJ/m2 (1500 Btu/ft2) prior to the initial thermal perfor- 
mance measurements. Series 2 is intended to determine the significance of the 
three-day preconditioning exposure specified in ASHRAE Standard 93-77. The complete 
set of baseline data is tabulated in Appendix C. All the measured data required to 
calculate and generate the efficiency curves are recorded for each data point. In 
addition, the time, wind, diffuse fraction, collector orientation and tilt are listed. 


Baseline data for all collectors measured are listed in Table 3 by intercept 
[optical efficiency, Fr(ta«),] and slope (heat loss coefficient, FRU). The curve 
fits and resultant values reported by the test sites were verified using the data 
points supplied by each test facility. 


The baseline data for the incident angle modifier consists of measurements made on 

most of the collectors at Sites 1, 2 and 3. Measurements were made on Collectors C 

and G along each of two orthogonal axes, i.e. parallel and perpendicular to the 
fluid-flow axis of the collector, because of special cover configurations causing 
non-symmetrical incident angle effects. The other collectors were not measured in 

the parallel direction. The test procedure describes two methods for incident angle 
modifier measurement. Test facilities No. 1 and No. 2 used a tracking mount and 
maintained a constant incident angle for each data point. Test facility No. 3 held 

the collector at a fixed tilt and orientation, and measured the efficiency when the sun 
was at the desired incident angle. The incident angle modifier can be expressed as: 


iL 
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where By = slope of the curve connecting data points plotted for Ky; versus 
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For the flat-plate collectors evaluated in this program, the incident angle modifier 


- 1). The baseline 
cos@ 


values of By for all measured collectors are shown in Table 4. Because of the 


sparsity of suitable clear days, no incident angle measurements were made at test 
facility No. 4. 


was essentially a linear function when plotted versus ( 
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DATA ANALYSIS 
5.1 General 


The data analysis was performed using the concepts and terminology recommended by 
ASTM Standard Recommended Practice E177 [12], where appropriate. A brief description 
of frequency, distribution, normal distribution, and the role of standard deviation 
is given in Appendix A. 


The test sites supplied linear and polynomial equations representing first and 
second order curve fits to the data points. The linear equations were used for the 
analysis used in this report to expedite the calculations and simplify the presenta- 
tion of heat loss coefficient data. Technical justification for the use of linear 
coefficients is provided in Section A-2, Appendix A. 


5.2 Thermal Efficiency 


The intercept and slope values listed in Table 3 are subject to two sources of 
variability: the random variability between individual collectors of the same type 
measured at the same site and the additional variability resulting from different 
sites. Each source is characterized by a standard deviation. The former is called 
the "within" (sites) standard deviation (s;). The standard deviation encompassing 
both "within" sites variability and variability among sites is referred to as 
"between'' sites standard deviation (sp). Detailed descriptions of the statistical 
analysis methods are given in References 1S A, gia! 15}, 


It is also possible to calculate the uncertainty, determined as a standard deviation, 
of the grand average thermal performance for each type of collector. Such a standard 
deviation, pertaining to an average, is generally called a standard error. Tables 5 
and 6 list the mean, standard error and the standard deviation "within" and "between" 
sites for intercepts and slopes, respectively. The standard deviations are listed 
both in their original form and as relative percentage values of the overall mean. 


There appears to be no relation between the standard deviations for any particular 
type of collector and the corresponding mean values. The variability exhibited by 
the standard deviations in a particular column of Table 5 and 6 is within the 
expected sampling errors of standard deviations based on the number of samples 
involved. 


As might be expected, the "between" standard deviations that include site-to-site 
variability are somewhat larger than those measured for within-site variability. For 
intercepts, the coefficient ofvariation (standard deviation expressed as relative 
percent of the mean value) is about 2.1 percent for "within", cvy,and 2.4 percent for 


"between", Cvp. The corresponding values for the slope are approximately 5.9 percent 
and 8.4 percent. 


Graphical depiction of the spread of individual data points from each test site and 
a comparison with the +4 percent to +5 percent random error band associated by Jenkins 
and Hill [7] with the instrumentation and represented by the dashed lines is shown 
in Figures 4a and 4b for Collectors D and B, respectively. It should be noted that 
no data were rejected on the basis of being an "outlier." The data points include 
variations in factors such as the measurement and calculation of collector area, 
flow rate, tilt angle, and orientation (tracking or fixed). 


7 


EFFICIENCY 


COLLECTOR B 
RESIDUAL STD. DEV. = 0.0202 


Oo1 02 03 04 05 06 07 08 
tf i-ta °C-m?2 
ew 


Figure 4a. Comparison of data point spread for double covered selective 
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Figure 4b. 


absorber type collector with uncertainty band 
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Comparison of data point spread for double covered nonselective 
absorber type collector with uncertainty band 
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5.3 Incident Angle Modifier (IAM) 


Although the total sample size for incident angle modifier (IAM) measurements is 
only one-half that for thermal efficiency, these data are important because of the 
lack of comparative IAM data and the interest in values for calculating all-day 
performance. Analysis of the data listed in Table 4 was performed for two individual 
test sites and combined for three test sites as shown in Table 7. The mean values 
for Collectors A and G exhibited relatively large differences for "within" data. 
The corrugated outer cover on Collector G probably influences the angular trans- 
mittance, and this combined with the difference in diffuse solar fraction between 
test sites (about 10 percent at Site 1 and about 25 percent at Site 2) could con- 
tribute to the large differences found for this collector. There is no apparent 
reason for the large difference for Collector A. The "within" test site standard 
deviations do not exhibit a consistent trend with collector type. 


Analysis of the IAM data to obtain "between" site standard deviation shows a large 
variation with collector type but the average coefficient of variation (expressed 
in percent) of 34 percent remains about the same with or without the extreme high 
or low values shown for collectors A, C, E, and F. Variations in IAM of this 
Magnitude cannot be attrivuted to product variability for the flat-plate collectors 
used in this program. 
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IMPACT OF DATA UNCERTAINTY 
6.1 General 


The uncertainty associated with commercial collector performance data is seldom 
published by the manufacturer and is not required to be reported by the test method. 
Three applications for typical thermal performance data have been selected for 
illustration of the impact of the data uncertainty discussed above, as follows: 


o Evaluation of collector degradation 
o Collector rating 
o System performance prediction. 


6.2 Collector Degradation 


One method of quantitatively evaluating changes in the material properties or thermal 
characteristics of a collector is to perform thermal efficiency tests before and 
after non-operational exposure testing. Depending upon the data uncertainty, thermal 
performance measurements could reveal incremental changes in the cover solar trans- 
mittance, plate solar absorptance or emittance,and changes in insulation thermal 
conductivity. Analytical studies of probable error resulting from instrumentation 
only [16] indicated that an uncertainty band of about 4 percent in the thermal perfor- 
mance could be obtained. This uncertainty in knowing the actual thermal performance 
corresponds to changes in the solar transmittance, absorptance or emittance of about 
0.10 and a change in back insulation thermal conductivity of about 50 percent for the 
collectors studied in this program. Consequently, changes in material properties of 
this magnitude would be required to clearly discern material property changes from 
the measurement uncertainty band for many flat-plate collectors. 


The influence of changes in the solar transmittance (t) of the cover material for 
Collector A is illustrated in Figure 5. The dark curves represent the calculated 
first order efficiency curve for t values of 0.90 and 0.82. The light curves repre- 
sent the instrumentation uncertainty associated with this collector and indicates 
that a change of 0.08 in t would be required to positively conclude that a material 
property change had occurred. An uncertainty of 0.025 at a nominal intercept value 
of 0.60 represents a 4 percent change in relative efficiency which is approximately 
equal to the 29 "within'' site coefficient of variation (4.1 percent) determined for 
intercept measurements reported in this study. Therefore, thermal performance 
measurements may not be sensitive to material degradation changes incurred during 
short-term (30-day) durability tests. 


6.3 Collector Rating 


Rating a solar collector on the basis of thermal performance is useful for several 
purposes. Among the primary ones of interest to the consumer are: 


o To rank collectors on the basis of all-day output 
o To determine the relative difference (%) in energy output for a 
specific application. 


One proposed rating method [17] for ranking collectors uses the ASHRAE Standard 93-77 
curve (efficiency and incident angle modifier) to calculate the collector daily 
thermal output for days with specified inlet temperature, solar radiation and 

ambient temperature. 


The influence of data uncertainty on all-day collector energy output was calculated 
for four flat-plate collectors with the properties and data uncertainties shown in 
Table 8. Calculations were made with a clear day solar radiation profile occurring 
at 40°N latitude, 0°C ambient temperature and a 60°C (140 F) collector inlet 
temperature. The data uncertainties used for the intercept, slope and incident 
angle modifier slope were determined from the baseline data analysis. Computations 
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were made for the nominal intercept and slope value, and the combination of uncer- 
tainties that would yield the lowest energy output, the greatest energy output, and 

a mixture of high and low uncertainties to give an intermediate energy output as 
shown in Table 9. Rankings based upon the outputs for each calculation series are 
shown not to change for the collectors evaluated. However, when comparing the energy 
output of a collector with a high combination of uncertainties to the energy output © 
of another collector with a low combination of uncertainties, a different ranking can 
be obtained. 


6.4 System Performance 


The prediction of a solar energy system's seasonal or annual performance with respect 
to meeting the heating or domestic hot water load can be calculated using system 
simulation programs. The f-chart method [18] is based upon a correlation technique 
derived from detailed modeling of a typical solar heating and hot water system in 

a number of U.S. cities with different combinations of component characteristics, 
sizes and loads. The collector performance, characterized by the first-order 
efficiency curve and the IAM slope, is used in conjunction with average local weather 
conditions and the site dependent load to determine the contribution to the load 
(fraction of the load) made by solar energy. F-chart calculations were made for a 
152 m2 residence occupied by a family of four and located in Madison, Wisconsin and in 
Washington, D.C. (see Appendix D for system description). The nominal collector area 
with a double cover selective absorber and a single cover flat black absorber, to 
provide 0.50 and 0.75 of the annual heating requirements for each location, were 
calculated as shown in Tablel0. The solar fractions for the same areas were also 
calculated using the high and low values of collector uncertainty obtained from the 
baseline analysis shown in Table 7. Differences of about +0.03 in solar fraction 

are attributed to the efficiency uncertainty for Collector 4 and about +0.045 for 
Collector 1. The 0.50 system fraction is somewhat more sensitive to the uncertainty 
values than the 0.75 system fraction for Collector 1. 


The influence of a +20 percent uncertainty in the IAM for either flat-plate collectoris 
shown in Table 10 to have less than a +0.01 effect on the annual system solar fraction. 
A sensitivity study of the relationship between the value of Bo and the normalized 
solar fraction for the solar heating and domestic hot water system example is shown 

for Collectors 1 and 4 in Figure 6. The solid curves are for the month of December 

and the dashed curves are for the annual performance in Madison, Wisconsin. Changes 

in the incident angle modifier slope of +0.05 or about +28 percent for a nominal 

value of Bo = 0.175 would result in a relative percentage change in annual solar 
fraction of about +3.0 percent for Collector 4 and +2.5 percent for Collector 1. 
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CONCLUSIONS ~ 


Statistical analysis of the thermal performance data sets obtained on eight liquid- 
heating, flat-plate solar collectors at four test sites broadly dispersed in the 
United States have confirmed that the total experimental uncertainty is about the 
same as the probable estimate of random error predicted for the use of ASHRAE 
Standard 93-77. Measurement error is believed to be the major contributor to the 
"within" test site variability (repeatability), and environmental effects a signifi- 
cant factor influencing the slope (loss coefficient) data obtained "between" 

test sites (reproducibility). No evidence of product variability was discerned for 
the various collector types. 


On the basis of first-order curve fits of the data points, an average "within" test 
site coefficient of variation, %¢vr, of 2.1 percent and "between" test site 
coefficient of variation, %¢vp, of 2.4 percent was calculated for the intercept 


[Fp (10) ]- 


Similar data for the overall heat loss coefficient, FpU;, indicated "within" (%cevy) 


coefficient of variation value.of 5.93 percent and "between" (%tvp) of 8.37 percent. 


The coefficient of variation (%cvp) for incident angle modifier data of 34 percent 
was shown based upon measurements at three test sites. Application of these uncer- 
tainty values to three application areas where thermal performance data have 
significant impact were shown as follows: 


o Material degradation equivalent to about a 0.10 change in cover 
transmittance or absorber solar absorptance and emittance 

o Variations in collector all-day thermal output of from +17 to 
+68 percent for specific winter operating conditions resulting 
‘in significant ranking changes 

o Variations in residential annual heating and domestic hot water 
solar fraction of +6 to 7 percent. 
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RECOMMENDATIONS 


The establishment of data uncertainty and measurement precision for commercial 
products requires careful evaluation of test procedures and continued vigilance of 
experimental facilities. It is recommended that a systematic evaluation of the test 
procedure be conducted within at least three test sites with replicate measurements 
(using the same "standard" collector types at each site). Collectors with stable but 
different material configuration designs should be selected. Some of the factors that 
should be considered are: 


o Absorbers with serpentine flow, headers with parallel passage and 
fully wetted 

o Curved covers vs. flat 

Selective vs. non-selective coatings 

o Relative high loss vs. low loss design (edge and back). 


eo} 


These measurements should be conducted with test apparatus that has been calibrated 
to the same field reference using a reference heat source or equivalent procedure. 

Measurements should be conducted under a range of local conditions as permitted by 

the test procedure (ambient temperature, wind, irradiance, diffuse fraction). 


Initial testing should be conducted on all collectors at one site using the same 
instrumentation. Series testing of groups of about four collectors utilizing 

one collector as a reference in all tests would minimize instrumentation uncertainties 
associated with the flow and solar irradiance measurements. The collectors would be 
measured once per month at each site over a six-month period (to include summer and 
winter) with the largest meteorological variations. At the conclusion of the program, 
the collectors could be returned to the site of the initial testing to corroborate 

the results if the uncertainties were larger than suspected from statistical principles 
of the sampling number. 
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APPENDIX A: DATA STATISTICAL ANALYSIS 
A-1 General 


The rectangles shown in Figure A-la constitute a histogram in which the base (A,B) of 
the rectangle ABCD represents a small interval of possible results of the measuring 
process. The area of the rectangle ABCD represents the relative frequency of values 
lying in the interval AB. As larger and larger data samples become available, the 
total range of values can be divided into smaller intervals. The histogram then 
generally approaches a smooth curve called a frequency distribution such as illus- 
trated in Figure A-1b. 


A precise measuring process is represented by a narrow frequency distribution. The 
smaller the precision, the wider the frequency distribution will be (see Figures 

A-le and A-1ld). The width of a distribution is measured by its "standard deviation," 
<= , which is a quantity expressed in the same units as the measured value. Thus, the 
standard deviation is an inverse measure of the precision. 


A type of frequency distribution that is of great theoretical and practical importance 

is the "normal" (or "Gaussian") distribution. In a normal distribution, an interval 
extending from a point located at one standard deviation below the mean to one located 

at one standard deviation above the mean, will correspond to a relative frequency 
(probability) of about 2/3 (68.27 percent). This is a +l sigma interval. The 

probability corresponding to a +2 sigma interval is approximately 95 percent and that 
corresponding to a +3 sigma interval is about 99.7 percent, as illustrated in 

Figure A-lc. The probabilities are measured by the area bounded by the curve on top 

and by the lower and upper limits of the interval on both sides. These probabilities 
pertain to the normal curve and one not necessarily correct for other types of distributions. 


A measuring process generally contains more than one source of variation (imprecision). 
Each source will have a frequency distribution with a definite sigma value. The joint 
effect of the various sources of error is a frequency distribution resulting from the 
superposition of the various individual curves. The standard deviation corresponding 
to this overall curve is the square root of the sum of the squares of the individual 
standard deviation. The square of the standard deviation is called the "variance." 
The central value, or "mean," of a frequency distribution is not necessarily the "true 
value" of the quantity measured. The distance between the "true value" and the mean 
of the distribution is called the "bias" of the method. Often the true value is not 
known and must be replaced by a "nominal" or "reference" value. A method with a large 
bias is said to be inaccurate or of low accuracy (see Figure A-le). 


A-2 Residual Standard Deviation 


When a curve (or straight line) is fitted to a set of data points, the vertical 
distance between the observed points and the fitted curve can also be measured in 
terms of a standard deviation. The smaller this "standard deviation of fit' (also 
called standard deviation of residuals or residual standard deviation), the better 
the fit. 


A comparison of the standard deviation of the residuals for the curve fits for each 
data point for Collectors C and E are shown in Figures A-2a and A-2b. Collector C 
had the lowest residuals (least deviation from a linear fit) and Collector E had the 
greatest deviation. Relatively little difference is indicated between the 

first- and second-order curve fits; consequently, first order curves were used to 
represent the intercept and slope. Table A-1 lists the residual standard deviation 
for curve fits to each collector data set. 
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EFFICIENCY 


Figure A-2a. 


EFFICIENCY 


Figure A-2b. 


01 


COLLECTOR E 
Ist ORDER Sy = 0.0176 
2nd ORDER Sx = 0.0163 


02 03 04 05 06 O07 08 
tfi-ta °C-m? 


| W 


Comparison of first and second order curve fits for FRP covered 


nonselective absorber type collector 


COLLECTOR C 
Ist ORDER Sy = 0.0055 
2nd ORDER Sy = 0.0036 


02 03 04 05 06 407 .08 


tfi-tg °C-m? 
ie aw 


Comparison of first and second order curve fits for one cover, 
heat trap, nonselective absorber type collector 
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APPENDIX B: DETAILED COLLECTOR DESCRIPTION DATA 


Table B-1 lists the important dimensions and material properties required for thermal 
analytical modeling. Most information was obtained from the manufacturers' 
literature. When the literature was not available, dimensions were taken from one 
collector. Handbook property values for thermal conductivity, refractive index and 
extinction coefficient were used when not specified by the manufacturer. 


Table B-2 lists collector area measurements made for each collector at Test Site 1. 
Calculations of the standard deviation and “within" test site coefficient of variation 
are shown. 


Table B-3 lists collector area measurements reported by each facility. Calculation 
of the standard deviation and "between" test site coefficient of variation is shown. 


38 


0:7W/Mx 
C 
233ed moT} eutjUuedizes sey kemosessede *poplTem—-yoqTIs st aqetd 12q108qy7- *pepuog einsseid st ajeqtd EO 


(Iv) | (1993s) } ~ (Ty) (tv) | (12228) (Tv) (tv) | (12298) 
Sz0°0 09°0 ZG 09°T 9L°0 O°T 79°0 8S°0 ce SA UB b ARIA! 


ZO (0,°W/M PATFIDEFFA) 


N 
Z0°0 | coro 720° z0°0 70°0 £70°0 | vor S€0°0 (ep antes Gone OT: 


ee | ee gre (err (ex T0°T 96°€ (4934T) eunToa 
9° ZST pe Sree eon | WN 41 IT (wm) [-Supoeds 
Bot» | Seo ceo 6ss'0 | €18°0 | 1s80 sexo | €s6°0 | (@) 4aPTM eaT300z3a 
mee | Sire eis ae ean eee men [wee] (m) yq3uey moTg 
Oe || Pez IP Sas oe Sa | 09€ o9¢ | sy | AatATaonpuog TewreUL 
Ber | eee Ae ro || OG a] aro so) ane | (mm) ssouyoTYL 
gs90°0 | gsgoc | |seoo'o =| deat | s600°0 JUSTOTJJE0) uoTIUTIXY 


Fesg — 
ieee! 


eae oe 
G600°0 €00°0 JUSTIOTFJFION VOTIIUTIXY 


cS °T 7ST XOPpUL <ATIOPAZOY 
ce ge Oe ee ee 
UOTIEUSTSEg 1OIIST TOD 


uOoTIe[NsSUT yorg 


(=) 


asesseg PIntT™ 


(1993S PTE) 
278Td Teqiosqy 


Zaao) ToeUUT 


XOPpUT SATIOPIFOY 


(mm) ssouyoTUL 


co 


(um) y eoeds 


IaA0) 129NnO 


ae 


uotsueutg 10 Ajiedoig queuodwmo9 


SUOTSUSUT( JUeUTIIeg pue seTIsedoig TeTAVIeW AOIOETTOD “[-d@ PeTIeL 


39 


Table B-2. Individual Collector Area Measurements by Test Site No. 1 


Coeff. of 
Vans 1G) 


a. 


Average Coefficient of Variation, Zev; = .328% 
(%cevy = 0.212% without 1.14 outlier value) 


Table B-3- Collector Area Measurements for Each Test Site 


Coeff. of 
Wiebe, (7%) 


oOo FP SC fH eC oOo & 


A 
B 
C 
D 
E 
19 
G 
H 


Average Coefficient of Variation, %Cvp = 0.94% 
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APPENDIX C: THERMAL PERFORMANCE DATA 


The test data and measured efficiencies are tabulated in the ASHRAE Standard 93 
prescribed format for all efficiency tests. The dataare identified by collector type, 
series number, and laboratory number (test site). In view of the present accepted 
practice in this country for building technology, common U.S. units of measurement 
have been used for the data tables. 
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APPENDIX D: DESCRIPTION OF SOLAR SPACE HEATING AND DOMESTIC WATER HEATING SYSTEM 
SIMULATED USING F-CHART 


Collector 


o Area: Sized for each collector in each city to provide 50% and 75% of the 
seasonal space heating and hot water load 


o Orientation: Tilted equal to the latitude, 0°, azimuth (south facing) 


FRU, 


(Btu/°F-£t2-h) 


One glass cover 0.670 -1.080 
Black paint absorber 

Two glass cover 0.642 -0.56 
Selective absorber 


Collector Loop Heat Exchanger 


fo) Collector Type 


o None 

Main Storage Tank 

o Capacity: 15 Btu/°F-ft2 (1.85 gal/ft2) 
Space Heating Load Heat Exchanger 


° €L Cin 


UA = 2.0 


Space Heating Load 


o Building Effective UA: 12912 Btu/°F-day 
(8 Btu/degree day-ft2 for a residence having a 
floor area of 1614 ft2) 
o Internal Heat Generation: None (Building Effective UA given above is assumed 


to be a net value including the effect of internal 
heat generation) 


Domestic Hot Water System 
o Daily Draw: 70 gal/day 
o Cold Water Temperature: 55°F 


o Hot Water Temperature: 140°F 
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